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Within the Bogoliubov-de Gennes linearization theory of quantum or classical perturbations 
around a background solution to the one-dimensional nonlinear Schrodinger equation, we study the 
back-reaction of wave packet perturbations on a gray soliton background. From our recently pub- 
lished exact solutions, we determine that a wave packet effectively jumps ahead as it passes through 
a soliton, emerging with a wavelength-dependent forward translation in comparison to its motion 
in absence of the soliton. From this and from the full theory's exact momentum conservation, we 
deduce that post-Bogoliubov back-reaction must include a commensurate forward advance by the 
soliton itself. We quantify this effect with a simple theory, and confirm that it agrees with full nu- 
merical solution of the classical nonlinear Schrodinger equation. We briefly discuss the implications 
of this effect for quantum behavior of solitons in quasi-condensed dilute gases at finite temperature. 
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I. INTRODUCTION 

The linear evolution of perturbations around known 
exact solutions is one of the most basic tools of theoreti- 
cal physics. It is invoked constantly in almost every field. 
While the perturbative approach is conceptually obvious, 
this does not make it naive or trivial. It is typically well 
justified as the first term in a systematic expansion in 
a small parameter of the full theory, and it can be very 
non-trivial indeed, in cases where the background solu- 
tion being perturbed is itself non-trivial. Hawking radia- 
tion PUS], as derived for linearized quantum fields prop- 
agating in fixed spacetime metrics with event horizons, 
is perhaps the most exotic example of this very common 
general scenario. From the very fact that the perturba- 
tions can be so dramatically affected by the background 
configuration, however, one must expect on the basis of 
Newton's Third Law of Motion that the perturbations 
should in turn evoke changes in background itself. This 
is known as back-reaction. 

In quantitative and technical terms, back-reaction is 
simply one component of the continuing perturbative ex- 
pansion beyond leading order; but it is often qualitatively 
and conceptually much more significant than any small, 
higher-order corrections to the evolution of the rest of the 
excitation modes. The paradigm of perturbations around 
a fixed background is so important that it often deter- 
mines our entire conception of a physical system. The 
fact that the background is not fixed may be merely a 
higher order perturbation in a mathematical sense, but 
it may often induce a non-perturbatively large change in 
how we think about the system as a whole. 

Because of their extreme controllability, Bose-Einstein 
condensates (BEC) of atomic gases [UJB] provide an ex- 
cellent test system in which to study back-reaction as 
a fundamental problem. Quasi-one-dimensional conden- 
sates can serve especially well in this way, because general 
solutions to the linearized problem may be available for 
them in explicit closed form, even for non-trivial back- 



grounds. One class of such non-trivial backgrounds in 
one dimension are the gray soliton solutions to the non- 
linear Schrodinger equation (NLSE) with repulsive non- 
linearity [7], which is known in its role as the mean field 
theory for a dilute Bose gas as the Gross-Pitaevskii equa- 
tion (GPE) ■ We have recently presented a complete 
set of exact solutions to the linearization of the NLSE 
(known in the many-body context as the Bogoliubov-de 
Gennes equation) around any gray soliton background 
|10) . This provides a rigorous basis on which to examine 
the back-reaction effects of perturbations upon solitons. 

As a qualitative fact it is well recognized that solitons 
are indeed affected by perturbations, including quantum 
and thermal fluctuations. The integrability of the one- 
dimensional NLSE with no spatial inhomogeneity makes 
the interaction between solitons and other excitations 
rather special, but even in this limit it is not entirely 
eliminated. In cases where integrability is broken, for in- 
stance by a trapping potential, interactions between soli- 
tons and perturbations are of course more general. The 
interaction between solitons and general time-dependent 
long wavelength hydrodynamic excitations was derived 
in [llj according to a multiple scale treatment within 
the classical one-dimensional GPE with slowly varying 
trapping potentials. Dissipation of solitons via scatter- 
ing with excitations when three-dimensional effects break 
integrability were discussed in [IS], and more recently in 
|13j . In Q31 [TS] soliton evolution was analyzed within 
the so-called truncated Wigner approximation |16j . in 
which the initial condensate state is sampled according to 
quantum statistics, but its dynamics is treated classically. 
A different mean-field-based approach was used in [PT] . 
based on a stochastic and dissipative Gross-Pitaevskii 
equation (GPE) [18] derived as an approximate treat- 
ment for the full many-body dynamics. Still other recent 
work towards identifying the effect on solitons of their 
embedding in the larger dynamical system of the quasi- 
condensate has numerically solved a fully quantum dis- 
cretized version of the problem, computing correlation 
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functions from which some conclusions about soliton be- 
havior may be drawn [191 ED] ■ 

The contribution of this present paper is complemen- 
tary to all of the above previous work. On the one hand 
we treat only the integrable homogeneous problem and 
reach one perturbative order beyond the linear. On the 
other hand, however, our results are analytical and ex- 
plicit, and they are based on exact solutions to the lin- 
earized problem combined with the exact constraints of 
momentum and number conservation. 

In this paper we analytically study the leading or- 
der dynamical back-reaction of traveling small-amplitude 
perturbations on one-dimensional gray solitons. Based 
on the analytical solutions to the Bogoliubov-de Gennes 
equation (BdGE) [TU], we first show that in the pres- 
ence of a gray soliton elementary excitations experience 
a well-defined phase shift within Bogoliubov approxi- 
mation, while the soliton is unaffected. Wave packets 
thereby undergo a position shift when passing the soli- 
ton. By including post-Bogoliubov terms and applying 
momentum and number conservation, we find a corre- 
sponding position shift of the soliton, and derive an ana- 
lytical formula for this deterministic back-reaction effect. 
Comparison with numerical integration of the GPE ver- 
ifies our analytical results classically. Our linear results 
are just as valid quantum mechanically as they are clas- 
sically, and the quantum version of our post-Bogoliubov 
calculations introduces only a trivial operator ordering 
question. Since the conservation laws on which our con- 
clusions are based remain valid quantum mechanically, 
therefore, our analysis can be straightforwardly extended 
to the quantum problem. 

The paper is organized as follows. In Sec. [IT] the gray 
soliton as a class of solutions to the GPE and and the ele- 
mentary excitations around the gray-soliton background 
as solutions to the BdGE are reviewed, with special fo- 
cus on their asymptotic behavior away from the soliton. 
In Sec. |III| we construct excitation pulses of finite wave 
length in the gray-soliton background by superposition 
of the elementary excitations and find a displacement of 
the pulses when crossing the soliton. In Sec. |IV| we de- 
rive an analytical formula for the back-reaction of the 
soliton due to the excitations, which needs inclusion of 
post-Bogoliubov terms, and compare it with numerical 
time propagation under the GPE. In Sec. [V] we summa- 
rize our results and briefly discuss some implications of 
our results. 



II. GRAY-SOLITON BACKGROUND AND 
ELEMENTARY EXCITATIONS 

We begin by discussing gray solitons as background 
field configurations [11] . These are a family of mean 
field wave functions parametrized by the soliton position 
xq, constant soliton velocity /3, and fixed background ve- 



locity v and asymptotic density /i = c 2 : 

fo(x) = e-' m(x ~ Xo) {\p + KtanhK(a; - x )} | 

k = - fi 2 ■ 

The constant k must be real, and consequently the soli- 
ton speed /3 cannot exceed the speed of sound c = »/ju 
of the background field. For increasing j3 the density 
depression around the soliton position x becomes shal- 
lower until the gray soliton reduces to a uniform back- 
ground for j3 — > The gray soliton can be con- 
sidered a localized object, since the local density |i/)o| 2 
approaches the asymptotic value /i exponentially rapidly 
for — xqI 3> 2. The phase of ip® likewise rotates within 
the localized region k\x — xq\ < 2 between asymptotically 
constant values. With a compatible choice of the gas ve- 
locity v, therefore, the soliton ([I]) also represents a suit- 
ably periodic wave function for a finite ring-like system of 
perimeter 2L, up to negligible errors of order 0{e~ 2KL ) 

The soliton background (fTl) is an exact solution to the 
NLSE/GPE in 1 + 1 dimensions (one dimension of space 
and one of time) [5T] . For repulsive point-like interaction 
among the particles and uniform potential this reads 

ift#B,t)= (-\dl+i{p-v)d x + \*p\ 2 -p\*p, (2) 

expressed in appropriate dimcnsionless variables and in 
a frame moving with velocity f3 — v relative to the labo- 
ratory frame. We have defined 

v 2 

fl = n + vj3 - — , (3) 

to shorten many equations. Thus, ([I]) represents the gray 
soliton in a frame co-moving with the soliton. In the 
context of Bose gases the NLSE ^ is obtained from the 
Heisenberg equation of motion for the second-quantized 
particle destruction operator field in mean field approx- 
imation. Since the soliton position Xq is kept constant 
by our adoption of the co-moving frame, we hereafter set 
xq = without loss of generality. 

A. Small amplitude excitations within Bogoliubov 
approximation 

The next step is to consider small perturbations around 
the gray soliton background. Small amplitude excitations 
of any wave length in a Bose condensed gas may be de- 
scribed within the Bogoliubov linearized approximation. 
We assume that the field ip can be written as 

tp(x,t) = ip (x) +eipi(x,t) , (4) 

where e is a small perturbation parameter and the soliton 
background tpo and the excitation field ipi are assumed 
to be of the same order of magnitude. We do a pertur- 
bation expansion in e and obtain the BdGE for i\>\ as 



3 



the linearization of the NLSE around the gray soliton 
background: 

id t Vi - (-\% + up - v)d x + 2\m 2 - #Vi + ^1 , 

(5) 

in the frame co-moving with the soliton. A solution of 
the BdGE can be obtained in terms of a normal mode 
expansion with respect to complex variables a k {t) , a k lt)* , 
labelled with mode index k. From the results in [10] one 
finds that 



function fix) 



ipi(x,t) = e 



dk(u k (x)a k (t) + vl(x)a* k (t)) (6) 



is a solution to the BdGE ([5]), where the time evolution 
of the complex variables is harmonic, i.e., 

a k (t) = a k e- in *\ a%(t) = a* k ^ , (7) 

and the mode function Uk{x), Vk(x) are given by 



fl k 



kn sech kx — 2(3fl k 



(k 2 ± 2fl k ) - + ifctanh K-X 



(8) 



Here the upper (lower) signs of =p and ± apply for u k 
(v k ). The frequency Q k = ~(3k + \[W[i + fik 2 is the 
Bogoliubov frequency for gray-soliton background and 
the normalization constant N k is chosen such that the 
mode functions are normalized to Dirac (5-functions ac- 
cording to |22) . From the form of the exact solution one 
can see that excitations are not reflected by the soliton: 
it is absolutely transparent. The mode functions u kl v k 
constitute the non-zero frequency part of a complete set 
of functions. For our considerations in this paper we will 
not need to introduce explicitly the modes with zero fre- 
quency, which require a slightly different treatment |10j . 

In what follows we will construct wave packets by spec- 
ifying the amplitudes a k as a Gaussian envelope, but be- 
fore doing so, we will examine the asymptotic behavior 
of the mode functions (|8| to the left and right of the 
soliton, since this will offer a basic insight into the non- 
trivial motion of perturbation wave packets in the soliton 
background. As will then become apparent in Sec. |HI| 
this is in turn the origin of the soliton back-reaction dis- 
placement. 



B. Asymptotic behavior of the normal mode 
functions 

To understand the essential feature of the asymptotic 
behavior of our u k and v k , it is helpful to define the 
corresponding asymptotic limits of i/jq itself. Because this 
notion of asymptotic behavior away from the soliton will 
be basic to our paper we introduce the notation for any 



/>0) = Jim /O) 
/<0) = J im /(») 
and = lim f(x) 



for x > 
for x < 



2 =F2« 



(9) 



as a compact form to refer to either of the asymptotic 
cases alternatively. Referring to i/jo(x) as given by ([I]), 
we abbreviate i/'o^ as simply ip> and recognize that = 
e" ira (i/3±K). 

For x > the gray soliton solution ([!]) is equal to 
if> > up to errors 0(e~ 2K > x >), whereas for x < it co- 
incides with ij) < up to 0(e~ 2li \ x \). Both these wave 
functions ip > and ip < represent uniform, though not 
real, solutions to the NLSE |2]), with constant density. 
In particular both tp > and tp< have the same density 
IV^I = \/ k 2 + ft 2 = y/JI = c, and the same phase gra- 
dient factor e~ lvx ; they differ only by a constant phase. 
It is easy to show that the finite frequency Bogoliubov 
solutions |8]) for the soliton background reduce, in the 
same e T2KX — > sense, to Bogoliubov solutions for the 
corresponding uniform backgrounds ij)^ |10| - The inter- 
esting point is precisely how the full solutions interpolate 
between these two asymptotic regions. 

We will therefore compare our exact Bogoliubov solu- 
tions for the soliton background to the corresponding Bo- 
goliubov solutions for a uniform, soliton-free background 
with the same modulus and phase gradient as ip^ , namely 

■0 = k 2 + /3 2 e~ lvx . We introduce the functions u k ,v k , 



u k ± v k 



1 



k 2 /2 
2^ \ fl k + Pk 



(10) 



which, when multiplied by phase factors e lkx and e~ ivx , 
are exact Bogoliubov solutions for the uniform back- 
ground solution $ [2"5] , 

Inserting the above definitions and applying some al- 
gebra, one finds that the Bogoliubov modes Q can be 
re- written as follows: 



_ gjtffc sgn(x) 



x sgn(fc) 
1^1 



v k ip^ 



where the phase 9 k is given by 

lh _ 2\(kn 2 - pn k ) 



Kk 2 



-<D{e- 2K W) , 

(11) 

(12) 



2i{kK 2 - pn k ) - K k 2 ' 

What we have thereby shown is that exact solutions in 
the soliton background reduce to uniform solutions on ei- 
ther side of the soliton, but with a particular fc-dependent 
phase shift between the two asymptotic solutions, across 
the intervening soliton. We observe, for instance, the 
limiting cases 



lim e Wk 

k— »±oo 

lim e i8k 
fc-s-0 



0> 
1 . 



(13) 
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This implies that for very large k (for which v k — > 0) there 
is no order e phase shift across the soliton in ip = ij)( S -\-eip\, 
while for small k this total phase shift is simply given by 
the background phase shift that is already present in ipQ. 
As a function of k, we can say that the excitation phase 
shift 9 k thus interpolates between two different senses of 
being trivial. When we consider the propagation of wave 
packets composed from our exact Bogoliubov excitations, 
however, we will see that k has non-trivial effects for all 
k. 



III. DISPLACEMENT OF EXCITATION PULSES 
A. Bogoliubov wave packets 

From the solutions ([6| we will construct a class of 
Gaussian wave packet excitations. In the limit of short 
wave lengths these will reduce to ordinary Schrbdinger 
wave packets of free particles with velocity large com- 
pared to the speed of zero sound in the condensate back- 
ground through which they move. In the limit of long 
wavelengths they will rather represent pulses of zero 
sound, broad enough to have well defined wavelength. 
For intermediate wavelength these packets are simply lo- 
calized travelling excitations of an intermediate nature. 

We define our wave packets by choosing 



O-k+Ak 



-Ak A 



2n 



(14) 



where the free parameters k and A define the mean in- 
verse wavelength and total spatial extent of the wave 
packet, respectively. For the following we define the di- 
mensionless wavenumber difference £ such that 



Afc = C/A 



(15) 



and assume that A is large enough that ( 14 ) is a Gaussian 
distribution sharply centered around Afc = 0. 

We explicitly compute the resulting excitation ([6| for 
k\x\ ^> 1, i.e., in the domain outside the soliton, by ex- 
panding u k+ ^/\ and i^, + £/ A , as they appear in the inte- 
grand of m powers of 1/A. Inserting Q into ([6} 
and expanding the integrand to first order in A -1 we find 
by integration 



<• i 



l^l A 



ei(fe _n fet± ^)(- fc + i *± 4 



+ e -^-n kt ±^)l- k _ i z±-, k 



A 2 

>± 
A 2 



z± = x - 



v k t± 



(16) 



(17) 



where we have introduced the fc-dcpcndcnt group velocity 

k 2 + 2c 2 



~dk 



sgn(fc) 



(18) 



and position shift 



At = - 



nk 2 



d6 k 



dk n k (n k + /3k) 



(19) 



such that the envelope of the excitation pulse in (16) 
depends on the argument z±(x,t) as defined. In (16) 



and in the following primes stand for derivatives with 
respect to the wavenumber k unless otherwise stated. 

We have hereby constructed a wave packet with well- 
defined wavenumber k and Gaussian envelope of breadth 
A, moving at group velocity v k . Since |/3| < c, for k > 
(k < 0) the group velocity v k is positive (negative) 
and the envelope of the wave packet moves in positive 
(negative) x-direction. In the |fe| 3> 1 regime v k — > and 
the excitation (16) reduces to a Schodinger wave packet. 



Eqs. (16) and (17) imply that the wave packet has 



the relative spatial displacement ±A k /2 between the two 
asymptotic situations at x ^ 0. A wave packet, which is 
initially located in the negative s-domain and has k > 0, 
travels in positive x-direction and eventually passes the 
soliton; afterwards it will be seen to have been pushed 
forward by an extra distance Afc in comparison to where 
it would have been if it had continued moving at con- 
stant group velocity, as it would have in the absence of 
the soliton. The wave packet displacement Afc is a mono- 
tonically decreasing function of \k\, and in the limiting 
cases obeys 



A = 



lim Afc = 

fc->±oo k 



c(c — sgn(fc)/3) 
4k 



(20) 
(21) 



The k —> limit in ( 20 ) can be understood in terms of the 



multiple scale theory of solitons interacting with general 
long- wavelength backgrounds [TT], though the analysis 
from that perspective is somewhat involved, and may 
indeed perhaps be most easily understood by reversing 
our present viewpoint, and considering the advancement 
of the long wavelength packet as back-reaction from the 
soliton's motion. The high k limit, however, can be easily 
understood semiclassically. 

When the wavenumber k is much greater than the in- 
verse healing length of the condensate (set to 1 in our 
units), one can neglect v k and the BdGE reduces to a 
time-independent Schrbdinger equation for u k |23j : 



(22) 



--di + i/3d x + V(x) ) u k (x)e lvx = Eu k ( X y 



In this equation E is the energy eigenvalue and V(x) — 
—2K 2 sedi 2 KX is the mean field potential 2|-0o| 2 exerted 
by the soliton-carrying condensate on short wavelength 
particles. In this high k limit we can solve (22) in the 
WKBJ semiclassical approximation and find 



u k e lvx oc exp 



ik 



2k 2 



dysech Ky 



(23) 
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while the energy eigenvalue is Ek = fc 2 /2 — k(3. The first 
summand in the exponent of ( 23 ) is equal to the result 



for constant background solution and the second sum- 
mand reproduces the phase shift (21) for large k in the 
limit x — > oo. This result (23) has the obvious classical 



interpretation that a particle accelerates and then decel- 
erates as it passes through the conservative potential well 
V(x), and thus emerges with its initial speed unchanged, 
but having taken less time to cross the well than it would 
have at constant speed. 



B. Numerical propagation of wave packets 

For comparison with the above analytical perturbation 
theory, we have also investigated the motion of the wave 
packets by numerically solving the NLSE pi), in cases 
both with and without a soliton present, and for vari- 
ous initial wave packet conditions. These results show 
that our analytical formula ( 19 ) for the soliton-induced 
wave packet advance Afe to leading order in perturbation 



amplitude e and wave packet momentum width A is 
accurate. 

In Fig. [T] we show a typical example for the propaga- 
tion of the wave packet excitation. We compare the re- 
sult for constant background ip to a soliton background 
[fi = v = 0.5 and = 1) with equal asymptotic 
density and gas velocity. The perturbation amplitude 
parameter is e — 0.02 <C 1. At initial time the cen- 
ter of a wave packet with mean wavenumber k = 0.7 
and spatial breadth A = 12 is placed at x = —80 in 
both backgrounds, such that at initial time the envelopes 
\ipO + eipi\ 2 are identical. Under time evolution the wave 
packets move in the positive ^-direction. The numerical 
integration was extended to the same final time for both 
backgrounds, at which point the excitation had in both 
cases passed well beyond the position where the soliton 
is in the case where it is present. We compare these nu- 
merically exact propagations of the wave packets to the 
analytically predicted motion of its Gaussian envelope, as 
given by the wave packet group velocity and the soliton- 
induced shift Afe. Our analytical envelope motion has 
also included order A~ 2 corrections to the group veloc- 
ity, as described in Appendix [AJ since our total evolution 
time is long enough that this small velocity correction 
produces a noticeable shift in the packet position at late 
times. This group velocity correction is the same whether 
or not the soliton is present, and all other order A~ 2 cor- 
rections are too small to be seen in our plots. 

We have also numerically confirmed the wave packet 
displacement Afe for many values of k and soliton speed 
j3 other than those represented in Fig. [T| Since for short 
wavelength packets, according to (21), the packet dis- 



placement soon becomes even shorter than the wave- 
length, it is clear that the soliton-induced displacement 
of short- wavelength pulses is an all but negligible effect. 
A soliton may be passed, however, by arbitrarily many 
excitation packets. The cumulative back-reaction effect 
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FIG. 1: (Color online) Wave packet displacement in a soliton 
background. In (a) the initial density of an excitation wave 
packet (dash-dotted line) characterized by k = 0.7, A = 12 
and s = 0.02 and center at x = —80 in a constant background 
(red lines) and in a soliton background (black lines) with the 
same asymptotic density and gas velocity is shown; the arrow 
indicates the propagation direction of the wave excitation. 
The soliton is located at x — and its other solution param- 
eters are /3 — v = 0.5 and yfp, — 1. The solid lines represent 
the density after the same propagation time for both back- 
grounds. In dashed lines the analytical envelopes including 
higher order corrections in A -1 before and after the propa- 
gation are depicted. Fig. (b) is a magnification of the gray 
square area in (a) and shows the final density of the constant 
background (red) and soliton background (black) with the 
corresponding analytical envelopes in dashed lines. 



on the soliton can thus become arbitrarily large. To this 
back-reaction effect we now turn. 



IV. SOLITON DISPLACEMENT AS A SECOND 
ORDER PROCESS 

At first order in e we have shown that the wave packet 
is pushed forward by the soliton, while the soliton itself 
is asymptotically unaffected. The reason for qualifying 
the soliton's invariance as 'asymptotic' is that for pertur- 
bations with wavelength longer than the soliton width, 
the terms in Uk(x) and Vk{x) that are proportional to 
sech 2 Ka; have exactly the effect of translating the soli- 
ton (since they are proportional the spatial derivative of 
ij>o(x)). The terms in Uk(x) and Vk{x) proportional to 
tanh kx likewise have effects on ip = ip + eipi that are 
indistinguishable, for small fc, from perturbations of f3 



G 



in ipo. By performing our time-dependent wave packets' 
Gaussian integrals, however, one can directly show that 
all these soliton perturbations vanish except when the 
wave packet envelope is near to the soliton. There is no 
lasting effect on the soliton, at order e, from Gaussian 
wave packets of the form we discuss. 

The soliton cannot possibly remain exactly unaffected 
by a wave packet, however, even asymptotically. A wave 
packet carries a finite total density perturbation, and 
thus makes a finite contribution to the motion of the 
system's center of mass. Momentum conservation and 
continuity (particle number conservation) together im- 
ply that the system's center of mass must travel at an 
exactly constant speed. The forward jump by the exci- 
tation pulse as it passes the soliton contributes a brief 
interval of acceleration to the center of mass. It is there- 
fore inconsistent with the exact conservation laws for the 
soliton to act on the wave packet in this way, without 
any corresponding back-reaction. 

This apparent center of mass jump is not a contra- 
diction to our exact order e results, because the wave 
packet only contributes to the center of mass at order e 2 . 
So what we have deduced from continuity and momen- 
tum conservation is a constraint on the order e 2 post- 
Bogoliubov corrections: they must cancel the jump of 
the center of mass that is contributed by the wave packet 
jump Afe. As we will see, the back-reaction displacement 
of the soliton is one of these demonstrably necessary e 2 - 
corrections. But to obtain the post-Bogoliubov soliton 
back-reaction correctly, we must consider the corrections 
at post-Bogoliubov order generally. 



A. Wave function at order e 

We expand the NLSE ([2| up to second order in e with 

tp = tpo{x) + etp 1 (x,t) + s 2 ip 2 {x,t) (24) 

where ip2 is assumed to be of the same order of magnitude 
as ipo and i>\- At order e 2 the NLSE is then given by 

id t </>2 = (~\dt + i(/3 - v)d x + 2|Vo| 2 -fi)ih + 4>l^2 

+ 2|</> 1 |% + ^ 2 . (25) 

It is of the same form as the BdGE ([5| except for the 
driving term on the right-hand side of Eq. (|25[) , provided 



by the soliton background and the first order contribu- 
tion tpi. The solution to the homogeneous part of (25) 
is therefore again of the form ^. To find a particular 



solution of ( 25 ) we expand the second order piece ij)2 of 
the condensate wave function in terms of a complete set 
of functions: 



i> 2 (x,t) = dl{ui(x)bi(t)+vf(x)^(t)} 



(26) 

+ q z {t)R z {x)+i Pz {t)S z {x) ^ 
where 6; (i) and fy* (t) are complex functions of time and 
q z {t) and p z {t) are real functions of time. The continuum 
mode functions alone {ui,v{\ do not constitute a com- 
plete set of functions; we must also include the discrete 
zero mode with negative mass m z — —4k 124] . associated 
with the spatial translation of the soliton, to achieve com- 
pleteness |10j . The discretely normalized mode functions 
in (26) are 







(iv + d x )ii 






_ idpipo/m z _ 



k 2 sech 2 kx 
(k + i/3(tanh nx + nx sech kx))/(4k 2 ^ 



(27) 



The representation of the discrete zero mode with real 
variables q z , p z , rather than the usual complex co- 
efficients related to quantum mechanical creation and de- 
struction operators, is necessary because harmonic oscil- 
lator raising and lowering operators (and their classical 
counterparts) are singular for the free-particle limit of a 
harmonic oscillator. 

Since we have explicit analytic solutions to the homo- 
geneous equation ([5]), it is straightforward to construct 
the Green's function to solve ( [25] ) . In fact we need sim- 
ply insert ([26]) into J25~l) 



and use the fact that ui and 
vi are orthonormal solutions to ^ to obtain the readily 
intcgrablc first order differential equation for the h(t): 



(iflt-n,)6i(*) 

= fdx[ut (2^1 1 



(28) 



A particular solution is always offered by 
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bi(t) 



2tti 



-e 2 



[(i'-if+(i"-if] 



+ 
+ 



. e -i(n,/+nj//)i 

e i(n,,+n,„)t 

e -i(n,/-n,//)« 
ft, - fv + £V 



da; [u*«i/(2^ uj" + "00^") + v* vy (2i/jqUih + ipovi")] 
Ax [u*vi>(2<ip u*„ + ij>o v i") + v * u H 2 ' l l ) o v l" + i>oU*»)] 
dx [u*ui>(2ip u*„ + ipov*,,) + u*v*,,(2ip vi> + tpo u l') 
v*u*,,(2tpoUv + ipoVi>) + v* l v v {2-4>lv* l „ + tpoup)]\ 



(29) 



+ 



(Here the prime on the wavenumber I does not indi- 
cate any differentiation, but simply distinguishes I' and 
I" from I as an integration variable.) Any solution 
bi(t) oc e _lSli * to the homogeneous equation can also be 
added to this particular solution, to satisfy initial condi- 
tions. Similar integrals provide exact particular solutions 
for q z (t) and p z {t). 



These integrals ( 29 1 can be evaluated exactly if need 



be, since all the x-integrals can be performed for any t; 
but by considering the e llx prefactors in itj and Vi, we 
can see that for early and late times, respectively before 
and after the wave packet has passed the soliton, the in- 
tegrals over I' and I" evaluate to wave packets in x that 
are Gaussianly concentrated well before or after the soli- 
ton. (This is obvious, inasmuch as we know that ?/>i is a 
wave packet that propagates through the soliton without 
reflection.) For t either before or after the packet passes 
the soliton, therefore, we can evaluate (29) to essentially 



perfect accuracy by replacing the ^{x), ui(x) and vi(x) 
functions in its integrals over x with their asymptotic 
forms for k\x\ 3> 1. This makes the x-integrals in (29) 



quite straightforward, for early and late t, allowing all co- 
efficients bi (t) , q z (t) and p z (t) to be computed quite sim- 
ply, except during the brief time while the wave packet 
overlaps with the soliton. (The behavior of these coeffi- 
cients during the overlap interval can also be determined, 
but the complicated expressions in this case do not seem 
to us to add any conceptual understanding, beyond the 
qualitative fact that all the co-efficients change during 
this time.) 

We can then obtain explicit expressions for the I inte- 
grals in tp2, for all t either before or after the wave packet 
has passed the soliton. In general these are again rather 
complicated, since ui and v\ are non-trivial; but it is again 
apparent from the integrand's form (for t before or after 
the soliton-packet overlap) that ip2 does not include any 
new wave packets, either reflected or otherwise emitted. 
Rather, ip2 includes just two kinds of contributions. 

Firstly, the single propagating wave packet is dressed: 
it acquires additional components, at order e 2 , that are 
either higher or lower harmonics of the order-e wave- 
The higher harmonic component is a co- 



traveling wave packet with mean wave number 2fc, while 
the lower harmonic component is a smooth pulse whose 
only spatial scale is the packet breadth A. Both of these 
harmonic packets travel with the primary packet, shar- 
ing its group velocity and dispersion; they do not sep- 
arate from it to move at the speeds at which isolated 
disturbances of their respective forms would propagate. 
Away from the soliton, this post-Bogoliubov dressing of 
the ipi wave packet reduces to exactly the same form that 
one finds if the calculation is repeated with a soliton-frcc 
background ipo, except that with the soliton the wave 
packet is displaced by the same as derived above: 



: 0F 



z|/2A 2 



e ~± 

^2 



-(Vk + (v + /3)Cfc) 



+ i^erf(^)((*/ + /?)% + c 2 a) 



v4 ast + o{\- 



where we have defined the quantities 

Vk = (\u k \ 2 + \v k \ 2 + u k Vk) , 
0c = {vku' k - u k v' k ) , 



(30) 



(31) 



and used the definition z± 



u k tT Afc/2 from (17) 



form 



In (30) we have explicitly displayed only the smooth 
dressing of the lower harmonic component; -0| ast denotes 
the higher harmonic terms proportional to the fast spa- 
tial oscillations e ±2lkx , whose explicit form will not be 
needed, for reasons that will be clear in our next Section. 
The lower harmonic dressing includes a phase pertur- 
bation with amplitude of order A -1 , whose spatial pro- 
file is the integral of the wave packet's Gaussian enve- 
lope. This phase perturbation thus has a profile pro- 
portional to the error function, which has the property 
that erf(oo) = erf(— oo) + 1. This means that a broad 
Bogoliubov-de Gennes wave packet at any finite k per- 
manently shifts the phase of tp behind it as it passes, by a 
constant of order e 2 /A. It is perhaps surprising that a lo- 
calized wave packet, even with high k, has in this sense a 
long range effect. But this is simply due to the fact that 



the sub-harmonic post-Bogoliubov components vary on 
the length scale A of the wave packet envelope, and so 
for A much longer than the healing length, these spa- 
tially smooth subharmonic corrections are in the hydro- 
dynamic regime of the nonlinear Schrodinger dynamics, 
in which density pulses are necessarily accompanied by 
phase steps. This subharmonic dressing effect is present, 
unchanged, in the absence of solitons. 

The second contribution in ip2 is one that persists in 
the vicinity of the soliton after the wave packet has passed 
it: the soliton is slightly displaced. It is otherwise exactly 
unchanged from its initial state, once the wave packet 
has passed it. This displacement is the back- reaction 
that we have been seeking, and in total it comes from 
two sources. Firstly there is a non-vanishing contribu- 
tion from the continuum of bi modes, which settles to 
a constant value even after the wave packet has passed, 



because the term in the integrand of ( 26 1 proportional to 



sech 2 ra has finite weight at small I. "Xnd secondly there 
is a contribution from the discrete soliton zero mode. 
The zero mode momentum coefficient p z (t) returns to its 
initial value after the wave packet has passed, and this 
initial value may be set to zero without loss of generality 
since an initial p z can be absorbed into the ft of ipo. But 
the soliton displacement amplitude q z (t) changes from 
zero to a non-zero value as the wave packet passes, and 
remains at this constant value thereafter. 

Both of these contributions to the soliton displacement 
back-reaction can in principle be computed directly, with 
q z and p z obtained by solving their own equations of 
motion, 

Pz = -2Re/dsi?:[2|7/; 1 | 2 Vo + V'i^] 

Pz 



4k 



+ 2Im/dzS:[2|?/> 1 | 2 Vo + VwS] ( 32 ) 



which together form the z-mode analog of ( 28 ) , namely 



the projection of (25) onto the zero mode subspace of 
tp2{x,t) that is spanned by R z {x) and S z (x). The mod- 
est difficulty of solving these coupled equations directly 
may be avoided, however, by using an indirect way of 
computing the soliton back-reaction, that is nonetheless 
just as exact as the direct approach, because it is based 
on exact conservation laws of the NLSE. 



B. Soliton displacement from conservation laws 

Having explicitly solved the NLSE in the domain 
k\x\ 3> 1 to second order in e 2 , we can now determine 
the back-reaction on the soliton from exact conservation 
of momentum and particle number. Since the NLSE con- 
serves both of these quantities exactly, its expansion in 
powers of e necessarily conserves them at each individual 
order in e. Both the wave packet dressing and soliton 
back-reaction contributions in i/> 2 contribute to total mo- 
mentum and particle number at order e 2 , and so from our 



exact expressions for and the dressing part of -02, we 
can use the conservation laws to infer the back-reaction 
contributions. 

To this end we define the center of mass Q and total 
linear momentum P of the system 



Q= dxx\^\ 2 , 



(33) 
(34) 



in the usual way. Using the NLSE ([2| (in particular its 
imaginary part, corresponding to particle continuity) and 
integration by parts, the exact equation of motion for the 
center of mass can be established as 



Q = P+[fi-v) / da#| 



(35) 



It also follows readily from the NLSE that the total 
linear momentum of the system P is exactly constant, as 



is the second term on the right hand side of ( 35 ) , which 
in the quantum gas context is associated with the total 
number of particles 



(36) 



N= dx\ip\* 



Consequently, the center of mass Q moves at exactly the 
same speed at all times. We can easily compute this 
speed to order e 2 during the evolution of our soliton and 
wave packet, by using the simple asymptotic forms for 
■01 and ip2 when the excitation pulse is well localized 
outside the soliton domain. We can simplify the expres- 
sions without loss of generality by choosing the reference 
frame in which the soliton is at rest (/3 = v), and thereby 
obtain simply 



Q = v k N + (D(e 3 ) , 



(37) 



where we define the two contributions Ni and N2 to the 
total number of atoms N = Ni + N 2 + 0(e 3 ) by 



N 1 = e 2 Jdx\i/jf\ 2 =e 2 



No. 



dx (ipotpf* + ipQipf 



(38) 
(39) 



2tt u\ + v\ + u k v k + iy k + (3)(v k u' k - u k v' k ) 



A 



(K + /3) 2 -c 2 )/(2c 2 ) 



In this case the velocity at which the center of mass prop- 
agates is determined by the group velocity v k of the exci- 
tation pulse and the integrated density of the excitation, 
i.e., the total number N of excited atoms it contains, 
including both Bogoliubov and (first) post-Bogoliubov 
contributions. (We can now see why the precise form 



last. 



of the higher harmonic dressing components e 



±i2kx 



i>2 

in (30) is not needed: it makes no contribution to N at 
order e 2 .) 
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Since the center of mass speed is exactly constant, for 
any times t± we must obviously have 



Q(t+) - Q(i_) = (t+ - t_)v k N + 0(e 3 ) . 



(40) 



But if we neglect the soliton back-reaction, and evalu- 
ate the left-hand side of ( 40 ) by inserting the asymptotic 
wave functions (16 1 and (30) in the definition (33), we 
find that it does not equal the right-hand side - if t_ de- 
notes a time before the wave packet interacts with soli- 
ton, and t + a time at which the packet has passed it. 
Explicitly, we find 



Q {>) (t + ) 



-n«) 



(*-) = (t+ - t-)u k N + A k N , (41) 



where with the superscripts > and < we indicate that we 
have computed the center of mass using only the corre- 
sponding asymptotic wave functions away from the soli- 
ton, which include only the packet and its dressing, but 
not the soliton back-reaction. As we have seen, the soli- 
ton has pushed the excitation packet forward by an extra 
distance and this has contributed an extra displace- 
ment of the center of mass at order e 2 by A k N . 

This cannot be the only contribution at order e 2 to 
the center of mass motion, for it is incompatible with the 
exact constant speed of the center of mass, as implied 
by the conservation laws. Hence we can infer that the 
soliton displacement contribution in ip% must supply the 
compensating correction, in order to preserve the exact 
result (40). Computing the left-hand side of (40) for the 
soliton background ([T]) , where at t + we account for a shift 
Ax in the soliton position by x — > x — Ax, we find that 
the soliton shift results in a center of mass displacement 
of 



Q (0) (i+)-Q tUj (i-) = -2kAx 



-o<°> 



(42) 



Here the superscript indicates that we have inserted 
the gray soliton solution tpo in (33). The factor —2k 



corresponds to the integrated linear density 'missing' in 
comparison to the otherwise constant background; in 
this sense it is proportional to the (negative) mass of 
the soliton. The center of mass displacement from the 
dressed excitation wave packet and the soliton back- 
reaction must cancel each other due to exact conservation 



of total linear momentum, and therefore Eqs. (40 1 to (42 ) 
yield 



Ax = A fc — = 



N _ fc 2 (A^! +N 2 ) 



2k 2tt k {tl k +j3k) 



(43) 



for the soliton displacement as back-reaction in post- 
Bogoliubov theory. Due to the negative mass —2k of 
the soliton, the soliton back-reaction displacement is in 
the same direction as the wave packet displacement: the 
wave packet drags the soliton with it a short distance. 
Furthermore one can see from the definition of the par- 
ticle numbers in ( 38 1 and ( 39 1 that the displacement is a 



Equation ( 43 ) is the main analytical result of this pa- 
per, and demonstrates the back-reaction on the soliton 
from wave excitations. It is determined by the finite shift 
of the excitations ( 19 ) in the gray soliton background, 



and the ratio of atom number atom in the excitation 
pulse and the soliton. The limit behavior of Afc, as dis- 
played in (21), implies that the back-reaction is maximal 
for k — and in the regime fc 3> 1 drops quadratically 
in k for fixed ratio N/2k. From the latter we conclude 
that short-wavelength particle- like excitations have a rel- 
atively small effect back-reaction effect on the soliton, 
whereas long-wavelength phase perturbations have the 
strongest interaction with the soliton. 



C. Numerical analysis 



The analytical formula (43) for the soliton displace- 



ment is now compared to numerical solutions of the 
NLSE for periodic boundary conditions using a split step 
method [25] . In this subsection we choose for the asymp- 
totic density and speed of sound \x = c 2 = 1, and place 
the initial soliton at x = 0. We emphasize that all re- 
sults presented above for the infinite system are valid for 
the finite system to order O (e _2fcL ); but a few compli- 
cations of detail are introduced by embedding the soliton 
and wave packet in a finite geometry. 

For periodic systems there must in general be a back- 
ground phase gradient ipQ oc e~' vx in the presence of the 
soliton, and the allowed values for the background flow 
speed v form a discrete set [TO] . Within this constraint, 
however, we can easily fine tune the system size 2L to 
achieve f3 = v, so that the soliton is then (initially and 
finally) at rest in the laboratory frame, and the soliton 
displacement Ax can easily be extracted from the density 
profile after time propagation. In the finite system, inte- 
grals over wavenumbers k naturally become summations 
over the discrete solutions to 



k cot kL 



Kk 



(44) 



according to [10]. By taking into account the freedom of 
an extra phase, which is allowed by the BdGE and the 
second order equation (25), we can place the center of 
the excitation pulse at initial time ti = at any position 
Xf, we choose a starting point well separated from the 
soliton. This is equivalent to replacing a k — > a k e~~ lkXi 



in (14). Furthermore we rescale our small perturbation 



parameter 



MK + v 2 ) 



(45) 



second order effect in wave packet amplitude. 



such that N\ is normalized to e 2 . 

A less obvious complication in the periodic realization 
of our problem is that, once the order-e 2 dressing is in- 
cluded, the total initial wave function that we have de- 
scribed in this paper heretofore is in general non-periodic, 
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FIG. 2: Normalized soliton shift Ai/e 2 as a function of the 
central wavenumber k for different soliton speeds j3. The data 
points represent averages of the normalized soliton shift over 
the pulse amplitudes e, which have been varied from 0.03 to 
0.3 in the simulations. The solid lines show the theoretical 
predictions with systematic inclusion of terms of order e 2 , 
while the dashed lines indicate the corresponding curves if 
second order contributions to the wave function tj) are not 
taken into account. The soliton speed was given by (a) /3 = 
-1/2, (b) /3 = -0.0058, (c) /3 = 1/3 and (d) /3 = 1/2. 



because of the phase step in (30). To realize on the ring 



an equivalent encounter between wave packet and soli- 
ton, therefore, we add to our initial state a BdG solution 
of order e 2 in the form of a hydrodynamic pulse [55] that 
cancels the total phase step around the ring, but then 
travels away from the soliton, and does not meet it or 
the main wave packet before the end of our numerical 
time evolution. 

We have solved the time evolution of the wave function 
under the NLSE for different soliton speeds /3 and for dif- 
ferent sets of the excitation parameters k and e, keeping 
A = 12 fixed. The parameters obey the conditions k > 
and kX 1, so that we have a wave packet of well de- 
fined wavelength moving in the positive direction. After 
the excitation has passed the soliton, the time evolution 
is stopped and the soliton position is extracted from the 
density profile of the final wave function [27 . For each 
value of the wavenumber k we have varied the amplitude 
parameter s between 0.03 and 0.3. The quantity Ai/e 2 is 
expected to be independent of e 2 from ( 38 1 , ( 39 ) and (43 1. 



Indeed, the relative deviation of this normalized soliton 
displacement from its average over the amplitudes never 
exceeds 0.2% for any k. Thus, the numerical solutions 
clearly confirm the soliton displacement to be an effect 
quadratic in the amplitude. 



In Fig. [2] the numerical results for Ax/e 2 , averaged 
over all chosen amplitudes e, are shown as circles, for 
different values of the wavenumber k and soliton speed 
f3. It is apparent that for fixed soliton speed the soliton 
displacement is a decreasing function of the wavenumber 
k, while for fixed wavenumber the displacement increases 
for increasing soliton speed. Furthermore Fig. [2] suggests 
that even excitations with central wavenumber k larger 
than the inverse healing length, which is one in our di- 
mensionless units, still have a noticeable impact on the 
soliton. The solid lines in Fig. [2] indicate the analyti- 
cal behavior, including all contribution to the density to 
order e 2 , for different soliton speeds (3. The analytical re- 
sult obtained from ( 43 1 and the numerical solutions show 



very good agreement, consistent with all our results being 
exact to the orders in e and 1/A that we have expressed. 
The qualitative behavior of the soliton displacement as a 
function of the soliton speed (3 and central wavenumber k 
seen in Fig. ^ can be understood in terms of the group 
velocity Uk of the excitation packet. It is obvious that the 
interaction time between the soliton and the excitation 
depends on the difference of the soliton speed and the 
group velocity. Hence, the soliton displacement increases 
with increasing soliton speed f3 and decreasing k. 

To demonstrate the effect of the wave packet dress- 
ing contribution in ip% on the soliton shift, we draw the 
analytical curve for the normalized soliton shift ignor- 
ing the contributions N2 ( [39] ) to the total number of ex- 
cited atoms N in Fig. [2j Without the systematic in- 
clusion of all second order pieces, the soliton shift is 
clearly underestimated, especially for low k < 1. This 
shows that the soliton back-reaction is in general a post- 
Bogoliubov effect, not only in the direct sense that the 
soliton motion zero mode is excited only at second or- 
der in e, but also in that the back-reaction is sensitive 
to the general post-Bogoliubov dressing of the other ex- 
citations. The number of excited atoms N 2 exactly con- 
tributes N/2 to the total number of excited atoms N for 
k = 0. The ratio N2/N is a monotonically decreasing 
function of k, and in the opposite regime fc 3> 1 one 
finds N 2 /N — 2^,/k 2 + O (fc~ 3 ); the excitation packet is 
dominated by the contribution N\. 



V. DISCUSSION 

In conclusion, we have demonstrated that a soliton un- 
dergoes a finite shift as back-reaction to passing small- 
amplitude wave excitations. Within linear Bogoliubov 
theory we have shown that excitation pulses are shifted 
forward by passing a soliton. Solving the NLSE beyond 
linear order, we have then derived the soliton displace- 
ment back-reaction from exact conservation laws. The 
size of the back-reaction effect can qualitatively be un- 
derstood in terms of the group velocity of the excitation 
and the soliton speed. The smaller the difference of the 
two velocities, the longer the interaction time between 
the excitation and the soliton, and the larger the back- 
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reaction. 

It is widely known that solitons are transparent to ex- 
citations, propagating through them without reflection; 
but our work here shows that this does not quite mean 
that solitons and excitations do not exchange momen- 
tum. What it means is that they have no net, perma- 
nent exchange of momentum. During the brief period 
when they overlap, however, they do exchange momen- 
tum temporarily. As an excitation passes through the 
soliton, it borrows some momentum from the soliton. It 
returns the full amount before leaving, but with the mo- 
mentum it briefly borrows, it travels some extra distance. 
The soliton correspondingly loses momentum briefly, and 
its position shifts slightly until it gains the momentum 
back as the excitation departs. Because of the solitons 
negative dynamical mass, this brief loss of momentum 
translates the soliton in the same direction as the ex- 
citation was advanced. The excitation is a trustworthy 
borrower of momentum, repaying its debt promptly and 
in full. The soliton is in turn a generous lender, since 
it charges no interest, but it is also a remarkably clever 
dealer, for after this peculiar transaction, both parties 
come out ahead. 

Our results shed some additional light on what ther- 
mal and quantum corrections to the simple GPE soliton 
must be like, because both thermal and quantum fluctu- 
ations can be represented as ensembles of wave packets 
like the ones we have discussed. On the one hand we 
can see that wave packets that are actually on top of 
the soliton tend to raise its temporarily. This means 
that if fluctuations really become strong, such that the 
soliton is continuously being passed by significant pertur- 
bations, then the darkest, narrowest solitons will indeed 
be smeared out, and become broader and grayer. On 
the other hand we can see that once a wave packet has 
passed the soliton, its lasting effect is only a translation. 
So in regimes with weaker fluctuations, we can expect to 
see something more like Brownian motion of the soliton, 
where it gets kicked back and forth at random intervals, 
while retaining its form. And we would expect the tran- 
sition between these two regimes to be gradual rather 
than abrupt, since both effects are always present, and 
which dominates depends only on how many perturba- 
tion packets are typically present. 



as well as the analog for Vk- The Gaussian integrals can 
then be written and solved in a compact way and with 
the definition of the phase <fi = kx ± 6^/2 — Q, k t we obtain 
outside the soliton range {n\x\ ^ 1): 



12 : 




exp 



2X1 



+ (u k O Vk, c.c.) 



(A2) 



The full expression of the excitation pulse at correspond- 
ing order is obtained by replacing Uk with Vk in the first 
line and taking the complex conjugate. With the second 
subindex i of ipu we indicate the order in A -1 we have 
kept. In (A2| we have defined \ 2 {jk) for an arbitrary 



function jk by 



*5 



A - — 



fk 
Jk 



i IKtT^ 



(A3) 



The effect of A 2 in comparison to ( 16 ) can be most easily 
seen if we assume e« 1, such that we can neglect terms 
of order e 2 in |-0o + £"012 1 2 ■ The envelope of the excitation 
is then determined by V ; o' ! /'i2+' ( /'o' ! /'i2- Since A 2 is complex 
the envelope of ^12 is corrected at second order in A~ 2 . 
Thus, the envelope is dispersive and the group velocity 
of the excitation pulse has increased to 



Appendix A: Group velocity to order A 

An elegant way of rewriting the integrand of ^ at 
second order is to use 



(u k + v k ) A 2 



+ o(\- 4 ) 



(A4) 



Uk+c/x = u k exp 



Xu k 2A 2 yu k \u k 

3~1 



xil + O 



It is important to notice that the change Svk = Vk — Vk in 
the group velocity, an effect of order 0(A -2 ), is present 
for both constant and soliton background. 
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